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ABSTRACT 

We investigate the age and metallicity distributions of bright globular clusters (GCs) in 
the candidate intermediate-age early-type galaxies NGC 3610, NGC 584 and NGC 3377 
using a combination of new Gemini/NIRI if '-band imaging and existing optical V,I 
photometry from Hubble Space Telescope data. The V—I versus I—K' colour-colour 
diagram is found to break the age-metallicity degeneracy present in optical colours 
and spectroscopy, as I—K' primarily measures a populations' metallicity. In addition, 
it is relatively insensitive to the effect of hot horizontal branch (HB) stars that are 
known to be present in massive old GCs. By interpolation between Simple Stellar 
Population model tracks we derive photometric cluster ages, metallicities and masses. 
In general we find that "metal-poor" GCs with [Z/H] < — 0.7dex are older than more 
metal-rich GCs. For the most massive GCs {M ^ 6 x IO^Mq) in NGC 3610 with 
available spectroscopic data, the photometric ages are older by ~ 2Gyr, and this 
difference is more pronounced for the metal-poor GCs. However, the photometric and 
spectroscopic metallicities are in good agreement. We suggest that this indicates the 
presence of a hot HB in these massive clusters, which renders spectroscopic ages from 
Balmer line strengths to be underestimated. To support this suggestion we show that 
all Galactic GCs with ^ 6 x lO^M© feature hot HBs, except 47Tuc. Using a recent 
observational relation between the luminosity of the most massive GC and the galaxy's 
SFR at a given age, we find that the galaxies' peak SFR was attained at the epoch 
of the formation of the oldest (metal-poor) GCs. The age and [Z /B] distributions of 
the metal-rich GCs are broad, indicating prolonged galaxy star formation histories. 
The peak ages of the metal-rich GCs in the sample galaxies are 3.7, 5.9, and 8.9 Gyr 
for NGC 3610, NGC 584, and NGC 3377 respectively. The peak value of the age and 
metallicity distributions of the GCs is correlated with the host galaxies' luminosity- 
weighted age and metallicity, respectively, indicating that the GCs can indeed be used 
as relevant proxies of the star formation histories of galaxies. 

Key words: galaxies: star clusters: general - galaxies: elliptical and lenticular, cD - 
galaxies: evolution - galaxies: formation - galaxies: individual (NGC 3610, NGC 584, 
NGC 3377) 



1 INTRODUCTION 

Globular clusters (GCs) are among the very few observable 
fossil records of the formation and evolution of galaxies. The 
nature of their integrated properties, closely approximated 
by a Simple Stellar Population (hereafter SSP), significantly 
simplifies the determination of their ages and metallicities 
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relative to that of the diffuse light of their parent galaxies. 
Hence they contain crucial information about the history 
of star formation and chemical enrichment of their parent 
galaxies (e.g., Whitmore et al.|1999 Goudfrooij et al.|2001 1 



A key development in this respect has been the discovery 
that the GC systems of many early-type galaxies have bi- 



modal colour distributions (e.g.. 


Kundu & Whitmore||2001 


|Larsen et al.||2001| |Peng et al.| 


2006|), providing clear ev- 



idence for the occurrence of a "second event" in the for- 
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mation of these systems. Such a bimodal distribution was 
actually predicted from merger scenarios of E galaxy for- 
mation by Schweizer (19871; Ashman & Zcpf (1992), who 
suggested that young GCs form out of chemically enriched 
gas during major gas-rich mergers, then redden and fade 
with time to produce the "red" population of GCs seen in 
old ellipticals. Since this prediction, populations of young 
and intermediate-age GCs have indeed been found in sev- 
eral gas-rich merger remnants (e.g. 



Holtzman et al. 



1992 



Schweizer et al.|[T996l [Miller et al.||1997| [Whitmore et al.| 



1999||Goudfrooij et al.||2001|200H 3; fCoudfrooij et al.|2007 l 
Although HST studies have provided much information 
on the bimodality of GC systems, the typically used V — I 
or g ~ z colours provide only limited constraints on the age 
and metallicity of GCs. This is due to the well-known age- 
metallicity degeneracy, where increasing age or metal con- 
tent have a similar effect in reddening the optical colours. 
The non linear relation between optical colour and metal- 
licity adds an additional complication in the interpretation 
of colour bimodality ( Yoon et al.||2006 |. In principle, spec- 
troscopy is a more accurate way to estimate age and metal- 
licity. However, in practice this is very time consuming even 
for an 8-lOm class telescope, and only 1-2 handfuls of GCs 
per galaxy have been studied in this way to date with suf- 



ficient precision (e.g.. 


Forbes et al. 


2001 


Goudfrooij et al. 


2001| IPuzia et al.||2005| |Woodley et al. 


20101 |Alves-Brito 



et al.|2011 1 



Another way to break the age- metallicity degeneracy is 
to combine optical and near-infrared colours. In the optical, 
hot evolved stars (e.g., horizontal branch (HB) stars and 
blue straggles) and young main sequence stars (if present) 
contribute significantly to the integrated light. They also en- 
hance the strength of high-order Balmer lines, which causes 
a bias toward younger ages when classically derived from 



such absorption lines (see more in Sect. 5.1 1. Conversely, for 
ages > 1 Gyr, the near-IR is completely dominated by light 
from red giant stars. Thus near-IR colours like I — K or 
J — K measure the typical temperature of the giant branch, 
which is known to be directly sensitive to metallicity with- 
out any significant age dependence (e.g., [Worthey et al.| 
1994). Relative to optical spectroscopy, the " optical-|- near- 



IR" imaging technique has the additional advantage that 
a larger number GCs can be detected in a reasonable ob- 
serving time per galaxy. This technique was pioneered for 



GCs in intermediate-age galaxies by Puzia et al. ( 2002 1 (see 
also [Goudfrooij et al.|200T| using ground-based (ESO/VLT) 
.ft'-band imaging in combination with Wide Field Planetary 
Camera 2 (WFPC2) V- and /-band data from HST to study 
NGC 3115 and NGC 4365, two nearby early-type galaxies. 



Puzia et al. ( 2002 1 found that NGC 4365 hosts a significant 
number of intermediate-age (2 — 6 Gyr old) GCs, whereas 
NGC 3115 only hosts " old" GCs. T he age and metalhcity 
values found for GCs by Puzia et al.| ( ,2002) were confirmed 
by subsequent spectroscopy of a handful of GCs in these 
two galaxies ( [Kuntschner et aL]|2002[ [Larsen et al.||2003| ). 
Finally, the complications with the non-linear relation be- 
tween optical colour and metallicity is safely avoided with 
optical-NIR colours (e.g.|Cantiello fc Blakeslee|2007||Kundu| 



fc Zepf|2007| [Spitler et al.|2008| . Hence, the optical-|-near-IR 
imaging technique is a valuable alternative to spectroscopy 
in studying GC age and metallicity distributions (see also 
Chies-Santos et al.|20lT |. 



The discovery of a significant number of intermediate- 
age GCs in a fairly "normal" elliptical galaxies like NGC 
4365 (Puz ia et al.||2002| |Kundu et al.| |2005', but se e [Brodie 
et al. 2005l|Chies-Santos et al.|2011l) arid NGC 5813 (IHempel 



et al.||2007^ , [2007j3) highlights the power of GCs in prob- 
ing major star formation events in galaxies. While such an 
intermediate-age population in NGC 4365 might not seem 
unexpected given the presence of a kinematically decoupled 
core, it has not been identified from the galaxy's integrated 
light. This raises the important question whether or not 
early-type galaxies, which show signs of past interactions 
or mergers, generally host intermediate-age, metal-rich sub- 
populations. 

In this paper we use the optical/near-IR technique men- 
tioned above to constrain the presence of intermediate-age 
population(s) of GCs in a small sample of nearby early-type 
galaxies that show signs of having undergone a dissipative 
merger several Gyr ago. Our target galaxies are NGC 3610, 
NGC 3377 and NGC 584 described further in the next sec- 
tion. 



2 THE TARGET GALAXIES 

The target galaxies were selected to have the following prop- 
erties at the time of proposal submission: (i) The galaxy is an 
E or E/SO galaxy that is nearby enough (whci < 2000 km s~^) 
to detect several tens of GCs in the K band in reasonable ex- 
posure times, (a) the galaxy has a rich GC system with high- 
quality optical GC colours measured with HST {V—I 



Kundu 



fc Whitmore|200T I whose distribution does not show a clear 
bimodality (the absence of such bimodality is expected in 
the presence of a significant subpopulation of GCs that is 
metal rich and of intermediate-age; see also |Kissler-Patig| 
et al.|1998 l, and (lii) measurement of (luminosity- weighted) 
age from integrated-light spectroscopy or photometry is in 
the range 2-5 Gyr. Table [l] lists relevant basic properties 
of the galaxies in our sample. Population properties of the 
sample galaxies from recent literature is discussed in more 
detail in Sect.[7T|. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Observations 

We performed near-infrared K'-hand (2.12 /xm) imaging of 
NGC 3610, NGC 3377, and NGC 584 with the Near-Infra- 
Red Imager and Spectrometer (NIRI) instrument on Gem- 
ini North ( Hodapp et al.|[2003 l. Our program ID was GN- 
2004A-Q-17. We used the f/6 NIRI camera imaging mode, 
which provides a plate scale of 0.117 arcsec/pixel and a field 
of view of 2 X 2 arcmin'^, similar to the field of view of the 
WFPC2 and ACS cameras aboard HST. 

In order to reliably detect the presence of both "blue" 
and "red" GCs, with S/N~ 5 similar to that of the 
WFPC2 data, the targeted limiting magnitude was the typi- 
cal turnover magnitude of GC systems of "normal" galaxies: 
Mv = — 7.3 mag (e.g., Kundu fc Whitmore 2001 Larsen et al 



2001 1, corresponding to Mk = —9.8 for a typical colour of an 
"old" (age of 14 Gyr) SSP with [Fe/H] = -1.5 (e.g. |Bruzual| 



fc Charlot|[2003l ). The exposure times to reach the limiting 
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Table 1. Basic properties of the sample galaxies. 



Galaxy Classification Ay (m - M)o D My (V-7')o {I-K)o AgCspcc [Z/H]spec SFRpcak 

RSA RC3 [Mag] [Mag] [Mpc] [Mag] [Mag] [Mag] [Gyr] [dex] [Mg/yr] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



NGC 3610 E5 E5: 
NGC 584 SOi E4 
NGC 3377 E5-6 E6 



0.03 32.71 ± 0.08=^ 
0.14 31.52 ±0.20=^ 
0.11 30.17 ±0.16'' 



34.84 -21.99 1.14 
20.14 -21.31 1.14 
13.87 -20.01 1.11 



1.96 1.7±0.1<= 
1.91 2.6 ±0.3'^ 
1.81 3.7 ±0.9'* 



+0.76 ±0.16^ 234 
+0.46 ±0.03^ 101 
+0.19±0.05<' 75 



Note. — Column descriptions: (1): Galaxy name. (2) and (3): Classification in RSA | |Sandage &: Tammann |1981[ l and RC3 | |de| 
[Vaucouleurs et aI.||199T^ catalogs, respectively. (4): Foreground reddening in Ay according to | |Schlegel et al.|1998| . (5): Distance 
Modulus according to reference given (see below) and the corresponding distance in Mpc in (6) and galaxy absolute magnitude in (7). 
(8) and (9): {V—I)o and {I—K)o colours according to |Michard| ( |200"5[ . (10) and (11): Spectroscopic age and [Z/H] according to reference 
given (see below). ( 12 ): Galaxy pea k SFR calculated in this study ( see Sect.|7.2|l 



Tonry et al. (2001 1; ' Harris et al. I 2007 1; ' Cantiello et al. 1 2007 1; ' Trager et al. 1 2000 I; 



Howell 



(20051. 




Figure 1. _ft:'-band Gemini/NIRI images of NGC 3610, NGC 584 and NGC 3377 from left to right, respectively. The outer four isointcnsity 
contours on the top panels indicate surface brightnesses of 5, 10, 50 and 100 X a above the background. Bottom panels show the galaxies 
with median filtered images subtracted. Different symbol types indicate globular clusters bluer (circles) or redder (diamonds) than 
I — K = 1.4mag or younger than BGyr (tilted solid line). NGC 3610 GCs with spectroscopic measurements in the literature are labeled 
and shown with larger circles. 



mas; nitudes (cf. Table[2| were estimated according to the dis- 
tance moduli and were such that NGC 3377 was observed in 
a single night, while NGC 3610 and NGC 584 required obser- 
vations to be spread over two and four nights, respectively. 
Table|2] summarizes the observations and night conditions. 
In Figure[l] we show the spatial coverage of the NIRI ob- 



servations. The character of the spatial distribution of de- 
tected GCs reflects the overlap between the NIRI and the 
HST/ACS,WFPC2 fields of view. Globular clusters shown 
with different symbols reflect their properties as derived in 
Section[5l 



4 /. Y. Georgiev et al. 



Table 2. Log of Observations 



Target 


Nights 


FWHM 


N 


Total Exposure 




yyyy-mm-dd 


arcsec 


imaf 


;es seconds 


NGC 3610 


2004-02-10 
2004-02-11 


0.43" 
0.35" 


64 
55 


3840 
3300 




2004-07-28 


0.45" 


35 


2100 


NGC 584 


2004-07-29 


0.46" 


34 


2040 


2004-07-30 


0.53" 


31 


1860 




2004-07-31 


0.41" 


38 


2280 


NGC 3377 


2004-03-02 


0.51" 


28 


1680 



3.2 Data Reduction 

Image reduction was performed with tasks of the GEM- 
ini/nirJ^ packages within iRAlj^ Due to the "first frame'j^ 
problem of the NIRI detector, all first exposures (dark, flat 
and on-/off-target images) of each new sequence were ex- 
cluded before basic image processing. 

Normalized flats were constructed from images taken 
with the calibration unit shutter closed ("lamps off') sub- 
tracted from exposures taken with the shutter open ("lamps 
on") and using short darks to identify bad pixels. Flat field 
images with the IR lamp on and off allow separation of 
the instrumental thermal signature from the sensitivity re- 
sponse. 

The near-IR sky level and structure varies on time scales 
of a few minutes. To account for such variations the imaging 
sequence of four on- followed by four off-target exposures 
provides close enough in time blank fields to derive the sky 
flat and correct the on-target exposures. 

Some of the NIRI frames contain an electronic pattern 
visible as vertical striping with a period of eight columns 
and not always present in all quadrants. Before including 
the affected frames in our image processing list we used the 
stand-alone python routine, nirinois e.p^which almost per- 
fectly removed the striping. The same image reduction steps 
were applied to the standard star observations as well. We 
discarded few images which had trails or double-peaked stel- 
lar profiles reflecting the loss of the guide star resulting in 
improper telescope guiding. 

All reduced science and standard star exposures of the 
corresponding nights were registered to a common coor- 
dinate system with niri/imcoadd based on GEOMAP geo- 
metric solutions. The flnal science image of each galaxy is 
derived from average combined individual images scaled to 
that of the most photometric night for the case of NGC 3610 



^ http: / / www.gemini.edu / sciops / data-and-results / processing- 
software 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ First frame problem arises because the NIRI array is not con- 
tinuously reset when idle. Thus, each new exposure after idle has 
different dark current level, probably due to image persistence 
after saturating the array (NIRI has no shutter). 
^ http:/ /staff.gemini.edu/'~astephens/niri/patternnoise/ 



and NGC 584, which were observed over few nights. The to- 
tal number of combined individual exposures for NGC 3610, 
NGC 3377, and NGC 584 is N=119, N=28, and N=137, re- 
spectively. The result is one deep stacked i^'-band image 
per galaxy. 



4 DETECTION, PHOTOMETRY AND 
CALIBRATION 

4.1 Detection and photometry 

Globular clusters in our ground based images would appear 
as unresolved point sources due to the relatively large NIRI 
pixel scale at the distances to the galaxies. For the nearest 
galaxy in our sample NGC 3377 at a distance of 10.8 Mpc we 
have a projected spatial resolution of 6pc/pix. This is twice 
the typical half-light radius (~ 3pc) of a GC (e.g. Barmby 



et al. 2007 Masters et al. 20101. NGC 3610 and NGC 584 are 



at much greater distances of 34.8 and 20.1 Mpc, respectively. 
Thus, even for the most nearby galaxy NGC 3377, the GCs 
on the ground-based images will appear as point sources 
with a stellar PSF, which determined the type of detection 
and photometry of GCs. 

To deteclrl as many GCs as possible to match with 
existing deeper HST photometry for these galaxies (Goud- 



frooij et al. 2007 Kundu & Whitmore 2001 1, we first cre- 



ated images representing the integrated light component 
of the galaxy by smoothing the combined images with a 
circular median filter with a radius of 8 times the seeing 
FWHM of the input image. The choice of the median filter 
radius resulted from tests designed to render the residual 
effect of point sources in the images negligible. This was 
done by measuring the magnitude difference between the 
original and median-subtracted images of several isolated 
non-saturated point sources. The original combined image 
was then divided by the square root of the smoothed im- 
age, providing an image with uniform shot noise charac- 
teristics (e.g., Goudfrooij ct al. 20071. This facilitates the 



uniform detection of point sources at a detection thresh- 
old of 3(T above the background, using the daofind task. 
The resulting object coordinates were used as input for PSF 
photometry performed on images from which the smooth 
galaxy light was subtracted. PSF models with radius 3 x 
the point source FWHM were created using several non- 
saturated point sources. PSF-fitting magnitudes were scaled 
using aperture photometry with a radius of 4 pix for the PSF 
stars. To estimate aperture corrections from curve-of-growth 
analysis, photometry of the PSF stars was performed with 
aperture radii of r = 2, 4, 12, 17, 21, 31 pixels. We found very 
similar aperture corrections for the PSF magnitudes among 
the three sample galaxies (—0.361 ± 0.008 mag). 

4.2 Photometric Calibration 

Photometric calibration was performed using photomet- 
ric standards (FS6, FS23, FS 34, FS 103, FSlll, FS 126, 
FS 130, FS 131 and FS 134) observed during the same nights 
as the target galaxies. Their instrumental magnitudes were 



^ All reduction and analysis has been performed with IRAF pro- 
cedures. 
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Table 3. Calibration Coefficients for Equation |[T]| 





ZPk' 


cfi 


c2 


NGC 3610 


-23.66 ±0.03 


0.07 


-0.025 ±0.024 


NGC584 


-23.68 ±0.04 


0.07 


-0.027 ± 0.021 


NGC 3377 


-23.68 ±0.03 


0.07 


-0.026 



the photometric properties of aU GCs with K' photometry 
which are discussed in the following. 



measured with the same aperture as their measured K' cata- 
log magnitudes in Leggett et al. (2006 1. Thus, we derive con- 



sistent transformation solutions from instrumental to stan- 
dard photometric system for all objects in our photomet- 
ric lists. To allow for optical-to-near-IR colour terms in the 
photometric calibration of the K'-h&nA data, we used the 
V-band magnitudes of the photometric standard stars in Ta- 



ble 1 of Hawarden et al. ( 2001 1, who compiled measurements 
from|Sand age fc Katem| (|1982[); |Lasker et al. (.1988 ); .Mc~ 



I ° — I I — — ^ , |- 

Cook fc S ion ( 198 7[);|Carlsberg Meridian Catalogue| ( |1989[ ); 
Leggett ( jl992| ); jLandolt (1992 1. Hence, the derived photo- 
metric transformations will take into account the V — K 
colour of the object. 

We evaluated a least squares fit to the following photo- 
metric transformation equation: 



mK' = K' + ZPk' + cl X X^' + c2 x (V - K') 



(1) 



where mK' is the instrumental magnitude, ZPx' is the 
zero-point, cl and c2 are the airmass and colour term co- 
efficients. The zero-point and airmass coefficients were fit- 
ted around their typical value^ ZPx' = 23.68 mag, cl — 
0.059 mag/airmass. However, a best fit (low rms and coeffi- 
cients uncertainty) was achieved for fixed airmass coefficient 
at cl = 0.07mag/airmas^ while the zero-point and colour 
term coefficients were kept variable. This resulted in stable 
solutions for the zero-point and the V — K colour term co- 
efficients, close to their typical values (cf. Tablejsf. Due to 
the only one standard star observed during the NGC 3377 
observing run, we fit the ZP^' using a constant value for 
the colour term, which is the mean value obtained from 
the NGC 3610 and NGC 584 calibrations. Excluding both, 
cl and c2 from the fit results in a ZPj^i = —23.73 mag, i.e. 
within the uncertainty estimate of the ZP^i value with the 
airmass and colour coefficients included. All transformation 
coefficients are summarized in TableO 

To convert the GCs' instrumental to standard K' mag- 
nitudes, l/-band magnitudes of detected GCs in our images 
were adopted from the HST/ACS and W FPC2 optical V, I 
photometry in |Kundu fc Whitmore] (|2001[ ) for NGC 584 and 
NGC 3377 and | Goudfrooij et al.| ( |2007f ~for NGC 3610. To 
match the NIRI and HST photometric lists, we transformed 
the ACS and WFPC2 GC pixel coordinates to the NIRI 
image pixel coordinate system by deriving geometric trans- 
formation solutions with GEOMAP and applying them with 
GEOXYTRAN IRAF routines. The transformation from in- 
strumental to K' magnitudes of our sample GCs was per- 
formed by evaluating calibration equation ([T|, taking into 
account the V — K' colour of the cluster. In Table|4]we list 



^ 'http:/ /www.gemini.edu/?q=node/10104 

^ 'I'he mean value found for a typical variation of tfie water vapor 
amount on Manua Kea ( Tokunaga et al.|2002| 
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Table 5. Completeness limits for object detection from X'-band 
completeness test. My calculated using a mean V — K' = 2.5 mag 
colour. 



K', [mag] 

Figure 2. Completeness curves for object detection as a function 
of K' magnitude for the three studied galaxies. Dotted horizontal 
lines mark the 50% and 90% completeness limits shown in Tablejs] 
Smooth solid lines show the two-power-law fit to the completeness 
distributions. 



4.3 Completeness in K-band 

To determine the completeness of the K-band observations, 
we used the ADDSTAR task in IRAF to add artificial sources 
to the final stacked image of each galaxy. To avoid intro- 
ducing increased crowding to the images while retaining a 
high number of artificial stars per spatial and magnitude 
bin for reliable statistics, we added 100 stars per image and 
repeated this step 100 time^ i.e. in total IC* artificial stars 
per galaxy. The images were then run through the same ex- 
tended galaxy light subtraction, detection and photometry 
steps. Detected and input artificial star lists were matched 
to compute the completeness as a function of K' magni- 
tude. The completeness curves are shown in Figure[2] With 
smooth solid lines we show a two-power-law fit to the com- 
pleteness values at each 0.5 mag bin. These functions were 
later used to assess completeness levels as a function of age 



and metallicity (see Sect. 5.31 

In Table[5] we summarize the 50% and 90% completeness 
values for the three galaxies. The absolute l/-band magni- 
tude (My) was calculated using the galaxies' respective dis- 
tance moduli (see Sect. 3.1 1, and a.V — K' = 2.5 mag colour 
for the GCs, which is the mean of the observed range of 
2.0 V — K' < 3.0 mag. Thus, at the respective distance 
to these galaxies we observe, with 50% completeness, just 
at the GC luminosity function turnover magnitude (GCLF 
ToM) for NGC 3377, and about 0.5 and 1 mag brighter than 
the GCLF ToM for NGC 584 and NGC 3610, respectively 
This corresponds to a luminosities greater than a few lO^L©, 
meaning that our JT'-band GG detection samples the bright- 
est and most luminous GCs. We note, that the currently best 
estimate of the distance modulus of NGC 3610 was derived 



~ 3yrs after the observations (see Sect. 3.11, which were 
planned with a closer distance to NGC 3610. This led to a 
lower completeness than planned for NGC 3610. 



* The value of the seed in ADDSTAR was different for each frame 



Target 


K' 


Mv 


Lv 




[mag] 


[mag] 


[lO^Lv,©] 




50% completeness 


NGC 3610 


21.71 


-8.50 


2.1 


NGC 584 


20.92 


-8.12 


1.5 


NGC 3377 


20.08 


-7.59 


0.92 




90% completeness 


NGC 3610 


20.51 


-9.70 


6.4 


NGC 584 


19.66 


-9.36 


4.7 


NGC 3377 


18.95 


-8.72 


2.6 



5 VIK' COLOUR, AGE AND METALLICITY 
DISTRIBUTIONS OF GCS 

As already mentioned in the Introduction, the combination 
of optical and near-infrared imaging photometry is a very 
efficient method to access the distribution of ages and metal- 
licities for an entire globular cluster system of a galaxy by 
means of a comparison to predictions of stellar population 
synthesis models. The specific power of this " optical -l-NIR" 
method is twofold: 

(i) for a population of coeval stars older than > 1 Gyr, 
near-IR colour indices like I—K or J — K mainly sample the 
metallicity sensitive temperature distribution of stars on the 
Worthey''1994). In combination with an optical 



RGB (e.g. 

colour like V — I, which is sensitive to both metallicity and 
age, this allows to break down the age-metallicity degener- 
acy present in the optical colour indices. 

(ii) Relative to optical colours and spectral indices, near- 
IR colours are significantly less sensitive to the presence of 
evolved hot stars such as blue horizontal branch (HB) stars 
and blue stragglers. Blue HBs can be present in GCs older 
than ~ 8 Gyr (starting with metal-poor GCs; see e.g. Lee 
et al.|2000 Yi|2003) , and the relative number of blue versus 
red HB stars in GCs is known to depend mainly on age, 
metallicity, and GC mass (e.g., Carretta et al.|2010 Dotter 



et al. 2010 2011 1. Hence, the presence of blue HB stars 



in GCs can cause incorrect age determinations when using 



optical photometry or spectroscopy (see also Schiavon et al 



2004j [Cenarro et al.||2007[ [Percival fc Salaris||2011[ ), while 
NIR colours are largely insensitive to their presence. 

In the following, our age-metallicity analysis will be based 
onV — I versus I — K' diagrams. 

In Figurejs] we present the V—I vs. I—K' colour-colour 
distributions of all GCs detected in K' in our sample galax- 



ies. Those are compared to Bruzual & Chariot (20031 SSP 



models using a canonical Chabrier (20031 IMF, for several 
ages (1, 2, 5, and 14 Gyr) and metallicity values [Z/Zq = 
0.005, 0.02, 0.2, 0.4, 1.0, and 2.5, equivalent to [Fe/H] = 
-2.3, -1.7, -0.7, -0.4, 0.0, and +OAf\ Figure [s] shows 



^ We note that the quoted photometric [Fe/H] is in fact [Z/H], 
i.e. the total content of metals relative to hydrogen given in the 
SSP models. However, the difference between [Z/H] and [Fe/H] 
is 0.03 dex < A[Fe/H]< 0.06 dex for low to high metallicities, i.e. 
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3 7 11 15 0.5 1 1.5 2 



(I-K)o 



N 



(I-K)o 



Figure 3. Optical-NIR (V — /)o vs. (I — K)q colour-colour distributions of globular clusters (solid circles) in NGC3610, 584 and 3377 
from top to bottom, respectively. Encircled solid symbols mark GCs with spectroscopic measurements in NGC 3610 | |Strader et al.|2003[ 
|2004| . Small dots are GCs with colour error greater than 0.1 mag. colour indices are foreground reddening corrected. Middle panels 
show histograms of (V— /)o colour distributions of the GC population detected with HST (with solid line) and those with NIRI /f'— band 
magnitudes (dashed line). Smaller solid circles with dashed line error bars show GCs with colour uncertainty larger than 0.1 mag. Left 
and right panels show a comparison between cluster colours and expectations for ages and metallicities from Simple Stellar Population 
(SSP) synthesis models ( [Bruzual fc Charlot||2003l l. Different line-types indicate isochrones of different ages for all metallicities covered 
by the models (left panels) and isometallicity SSP tracks for ages > 0.5 Gyr (right panels). The arrows in the bottom panels show the 
direction of increasing age and metallicity, as well as a reddening vector for Ay = 0.5 mag. 



much smaller than the measurement error. To avoid confusion in 
the following we adopt the [Fc/H] annotation for the photometric 
metallicities as well. 



that the three galaxies generally show a large spread in age 
and metallicity among their brightest GCs. A notable sub- 
population of roughly coeval bright GCs can be also seen 
for NGC 3377 (see below). The I - K' colours of the ma- 
jority of the "blue" population of GCs as judged from the 
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optical data (i.e., V — I < 1.0) turn out to be consistent 
with old ages (> 5 Gyr) and low-metallicities on average 
(-2.3 < [Fe/H] < -0.7dex), while the /-A" colours of the 
optically redder clusters tend to indicate higher metallici- 
ties and younger ages ([Fe/H] > — 0.7dex and ~ 3 — 6 Gyr, 
respectively) . 

While the overall colour-colour distributions are in- 
dicative of a large spread in age and metallicity among 
the brightest clusters, there is a discernible population of 
roughly coeval GCs in NGC3377 (i.e., the GCs between the 
2 and 5 Gyr isochrones in the lower left panel of Fig.|3|. From 
our data alone, there is no indication for a spatial correlation 
among these clusters (cf FigjlJ . While such a subpopulation 
is not readily seen in the other two galaxies, this may be due 
at least in part to the photometric completeness limits (cf. 
Table|5|. The middle panel histograms in Figurejsjshow that 
our NIRI imaging detects 61% (64/106), 73% (81/115) and 
8% (50/611) of the NGC 3377, NGC 584 and NGC 3610 GCs 
detected with HST, respectively. Although the NIRI imaging 
samples 73% of the NGC 584 GCs observed with WFPC2, 
the relatively large photometric errors cause a larger scatter 
in their VIK' colour distributions, which prohibits a clear 
detection of any underlying subpopulation of coeval GCs 
similar to that observed in NGC 3377. 



5.1 Photometric ages and metallicities 

Employing the good resolution in age and metallicity pro- 
vided by the optical-NIR {V, /, K) colour indices, we derive 
photometric ages and metallicities for the GCs in our sam- 
ple. This is done by a standard minimization interpola- 



tion between SSP model tracks of Bruzual & Chariot ( 2003 1. 
A brief comparison with other SSP models and the reasoning 
why we chose [Bru zual fc C hariot] models for this purpose is 
presented in Section [5. 2[ The derived ages and metallicities 
are the weighted average of the age and metallicity of the 
two nearest tracks. The weights include the distance to the 
respective model and the photometric errors. This way we 
also calculate an upper and lower error to the photometric 
ages and metallicities. We list the derived properties of all 
GCs in our sample in Table|4] 

About a dozen of the brightest GCs in NGC 3610 have 
available Lick index analysis of optical spectroscopy ( |Strader| 



et al. 2003 20041. Strader et al. qualitatively discuss the 



clusters ages, but they do not explicitly provide ages for each 
cluster. Thus, using their Lick index strengths, we calculate 
ages for the GCs with the GONZO code ( Puzia et al.||2002 1, 



using the models of Thomas et al. (20031. Of those GCs, 



eight are within the NIRI field of view for which we can 
compare our ages and metallicities with their spectroscopic 
values. In Figure|4] we show these differences (in the sense 
spectroscopic minus photometric) as a function of I — K' 
colour. Errors in the top and middle panels in Figure |4] are 
the square root of the sum of the squares of the spectroscopic 
and photometric errors while errors in the bottom panel are 
their photometric uncertainties. 

The least squares fit to the data (dashed lines in Fig.|4| 
show that spectroscopic [Fe/H] and photometric metallici- 
ties [Z/H] derived from the VIK colours agree very well, 
keeping in mind that [Fe/H] - [^/H]= 0.03 to 0.06 dex (cf. 
Sect. [5] above). In contrast, the photometric age estimates 
of the most massive GCs in NGC 3610 are older than their 
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Figure 4. Differences between spectroscopic and photometric 
metallicities (top) and ages (middle) for NGC 3610 GCs. Dashed 
lines are the least squares fit. Bottom panel: photometric ages 
and metallicities for NGC 3610 GCs with spectra. The top panel 
shows that VIK derived metallicities are well consistent with their 
spectroscopic values. On average, photometric ages are older than 
spectroscopic ones by ^ 2 Gyr (middle panel), and the difference 
between them is more pronounced for low-metallicity GCs than 
those with higher Zs, possibly due to the presence of blue hori- 
zontal branch (see text in Sect. |5.l1 . 

spectroscopic values by about 2 Gyr on average. The age 
difference is larger for low-metallicity GCs than for higher 
metallicity ones. Such a difference can in principle be caused 
by the presence of blue horizontal branch (BHB), helium- 
core burning stars in these clusters. It is well known that 
spectroscopic luminosity-weighted ages, derived from the 
classical H^j vs Mg2 or Mgb absorption line indices, can 
be biased toward younger ages, especially for metal-poor 
GCs, because the high effective temperature of BHB stars 
{Tsf f > 9000 K) and blue stragglers enhances the strength of 



Xin et al. 



g- 


Lee et al.||2000 


Schiavon et al.j 


2004 


|Ocvirk||2010| |Percival & Salaris| 


2011 



2011 1. An example for such "age bias" are the 



metal-poor Galactic GCs NGC 3201 and NGC 5024 featur- 
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ing BHBs, whose spectroscopically derived age (from H/3 vs. 
Mg2) is ~ 8Gyr ( Perina et al.|[20TT |, while their ages from 



deep HST CMD fitting is > 12Gyr potter et ar.||2010 |. 
Furthermore, nearly all massive {M > 3 X 1 0= My < 
—8.9 mag) Galactic GC have hot HBs (e.g. Recio-Blanco 



et al. 2006 Lee et al. 20071, with the exception of the 



metal-rich GC 47Tuc. This may be due to its low central 
escape velocity (wcsc) for a massive cluster at that metal- 
licity, unhke the metal-rich NGC 6388 and NGC 6441 which 
do have BHBs and larger Vcsc (cf. Fig. 5 in Georgiev et al. 



2009 1 and thus better able to retain/accrete processed stellar 
ejecta. Puzia et al. ( 2002 1 performed a Lick index analysis 



for metal-rich Galactic GCs and they note that NGC 6388 
and NGC 6441 (which host BHB stars) show a stronger H/3 
index than the other metal-rich RHB GCs. However, |Schi-| 
avon et al. (20041 estimate ages of ~ 8 and ~ llGyr from 



Lick index spectroscopy for these two clusters, which indi- 
cates that spectroscopic age-dating of metal-rich GCs can 
be less biased by the HB morphology. This is what we do 
indeed observe in the middle panel of Figure|4] for W 6 and 
W31, two metal-rich GCs in NGC 3610: their photometric 
and spectroscopic ages are more consistent (see also Sect.[6|. 
More pronounced HB effect is expected for metal-poor GCs, 
which we indeed observe for the two metal-poor GCs in the 
middle panel of Figure]!] This indicates that in spite of un- 
certainties inherent to SSP models, the most likely cause 
for the mean difference between photometric and spectro- 
scopic ages for these GCs is the presence of hot HB stars, in 
particular for the metal-poor GCs. 

The validity of the method of deriving ages and metal- 
licities from V, I, and K photometry used in this paper is 
supported by the observation that dereddened V—I and V—K 
colours of old GCs in the Milky Way and M31 have been 
shown to be consistent with lBruzual fc Charlotl SSP model 
predictions for "old" ages (> lOGyr) throughout the range 
of [Z/H] sampled, (see Fig. 5 of [Puzia et ar][2002| . The ap- 



same M31 GCs data as in Puzia et al. 


(20021, i.e. V",/ magni- 


tudes from|Barmby et al.|(|2000||2001 


1 , augmented with new 


spectroscopic measurements from Beasley et al. ( 2005 1 and 



Caldwell et al. ( 2011 1, the latter providing updated E{B—V) 
values. Even though the E{B—V) are uncertain (a~0.1 mag, 
see Caldwell et al. 2011 1, the V— 1, 1~K colours are consistent 
with thelB ruzual Sz Charlot| SSP model predictions, even so 
for GCs massive (M > 6 x lO'^ Mq) and £(B-V) < 0.3 mag. 
We derive their photometric ages and metallicities, which 
compare well with their spectroscopic values from |Caldwell| 
|et al.| ( |2011[ l, while spectroscopic ages for the oldest GCs 
{Agev IK >8Gyr), on average, are younger, i.e. showing the 



BHB effect. Due to this problem, Caldwell et al. (2011 1 does 



not derive spectroscopic ages for the majority of the GCs 
and assume age of 14 Gyr. We consider only GCs with spec- 
troscopic ages in Figure[5] 

5.2 Comparison to various SSP models 

It is beyond the aims of this paper to test different SSP 
models and investigate the nature of the differences between 
them. However, in order to provide a qualitative overview 
of how different synthesis population models compare to 
the observed colours of NGC 3610 GCs with available spec- 
troscopy, we compare our data with the popular SSP models 
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Figure 5. Comparison between VIK photometric and spectro- 
scopic ages and metallicities for M31 GCs. Top; Brightest (most 
massive) and less reddened {E{B—V) < 0.3 mag) M31 GCs, as in- 
dicated in the legend, which wc consider having the most reliable 
photometric and spectroscopic measurements. VIK magnitudes 
arc from [Barmby et"~al. (2000 2001) and spectroscopic ages and 
metallicities from |Beasley et al.| ( j2005) l and |Caldwell et alT] | |2011| . 
Colours were dereddened using E{B—V) values from these stud- 
ies. In spite of the large photometric and reddening uncertain- 
ties, we derive similar photometric and spectroscopic metallici- 
ties (middle panel), while spectroscopic ages for the oldest GCs 
{AgeYiK > 8 Gyr) are younger on average (bottom panel), sig- 
naling the BHB effect. In both panels, dashed horizontal lines are 
the least-squares fits to the data, with coefficients shown in the 
legend. 



of BC03 ([Bruzual fc Charlot|2003[), Galev ( [Anders fc Fritze- 
Alvensleben||2003| [Anders et al.||2009|), Vazdekis/MILES 



(Vazdekis et al. 



al. 



2000 Rai- 



2005 I SSPs 



2010|, TSPoT ( [Brocato et 
mondo et al. 2005 1, and the Maraston (Maraston 
for blue and intermediate-to-red HBs. Results are shown in 
Figure[6] As seen, the different SSP models yield significantly 
different results in terms of ages and metallicities. This is due 
in part to the different treatment of key stellar evolutionary 
phases, e.g. the TP-AGB and post-AGB phases, HB stars, 
etc., which are relatively poorly understood and empirically 
constrained. In addition, the different SSP models use dif- 
ferent stellar atmosphere spectral libraries (theoretical or 
empirical), and different stellar evolution prescriptions. The 
effects of different choices of the IMF on the integrated VIK' 
colours are, however, smaller than the observational errors. 
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Figure 6. Comparison between V — I , I — K' colours of GCs in NGC3610 with colour errors < 0.1 mag and different SSP synthesis 
models (see text in Sect. |5.2l . With circles are GCs which have spectroscopically derived ages and metallicities. In the top left panel 
with arrows are indicated the direction of increasing age and metallicity. With labels are indicated the 1,2,5 and 14Gyr isochrones 
and isometallicity SSP tracks with labeled values. All models are with canonical IMF (Kroupa 2001, Chabrier 2003). MILES SSPs are 
with the revised Kroupa IMF, which accounts for effects from unresolved binaries. Only Br uzual fc Chariot SSP models provide K' 
magnitudes, however, the K — K' difference is from —0.02 to 0.007 mag for high to low metallicities, respectively, which is negligible 
to account for differences between data and models. The Maraston models include the blue and intermediate-red horizontal branches 
tracks, which differ for ages > 5 Gyrs. 



Overall, our VIK' data is best represented by the BC03 SSP 
model, especially for low metallicities and old ages. This is 
why we use BC03 to derive and discuss ages and metallicities 
for our clusters in the following. 



5.3 GC age and metallicity distributions 

In Figure[7|we show the VIK' photometric age and metal- 
licity distributions for the GCs in the three target galaxies. 
Employing the good separation in the I — K colour in the 
SSP models (cf e.g. Fig.|3| between metal-rich and metal- 
poor and older versus younger populations, we separate the 
GCs in two subpopulations, namely redder and bluer than 
I — K — 1.4 mag, corresponding to [Z/H] ~ — Idex for a 
population older than about 5 Gyr. Within we then es- 
timate the highest-probability age for the metal-rich GCs 
using a non-parametric probability density estimator with a 



R is a language and environment for statistical computing and 
graphics | jhttp://www.r-project.org| . 



Gaussian kernel with size of 2 Gyr. Those probability den- 
sity estimates are shown with solid lines in Figure[7]and the 
peak values for the age distribution of the metal-rich GCs 
(AgCiod) are shown with labels. It is seen that AgCrcd in- 
creases from ~ 4 Gyr for NGC 3610, to ~ 6 Gyr for NGC 584 
and ~ 9 Gyr for NGC 3377. To assess the statistical level of 
similarity between the age distributions of the three galax- 
ies, we consider that one can not expect Gaussian (or any 
other anticipated functional form) for the cluster age distri- 
butions. Formally, this precludes the use of the Kolmogorov- 
Smirnov test which is commonly used for comparing statis- 
tical similarity between two samples. We therefore use the 
Mann-Whitney U testp^ which yields that the distribution 
of GC ages in NGC 3610 is different from those of NGC 584 
and NGC 3377 with a probability of 75% and 98%, respec- 

1^ The Mann-Whitney U test (a.k.a. Mann- Whitney- Wilcoxon 
rank-sum test or u-test) is a non-parametric statistical hypothesis 
test which calculates a measure of the difference between two 
independent samples of observations. It is similar to the t-test, 
but is more robust because it compares the sums of ranks and 
thus is less sensitive to the presence of outliers. 
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Figure 7. VIK photometric age and metallicity distributions 
(left and right panels) for all GCs in the three post-merger rem- 
nant galaxies, i.e. NGC3610, NGC584 and NGC3377, from top 
to bottom respectively. Light solid (gray), solid (red) and dark 
solid (blue) histograms show all GCs, redder and bluer than 
I — K' = 1.4 mag, respectively. Solid lines show the probability 
density distribution for all GCs. The highest probability density 
age for the red, metal-rich GCs, shown with label for each galaxy, 
increases from top to bottom. Dashed lines indicate the number 
density distribution corrected for 50% incompleteness (details in 
text in Sect . |5. Open solid line histograms present GCs with 
mass larger than the mass at 50% incompleteness at 14 Gyrs (cf. 
also Fig.[9]|. 

tively. The age distribution of GCs in NGC 584 diflFers from 
that in NGC 3377 with 98% probability. Thus, the brightest 
GCs in NGC 3610, NGC 584 and NGC 3377 have statisti- 
cally significant different age distributions from one another 
and the peak age of the metal-rich subpopulation increases 
from NGC 3610 to NGC 584 and NGC 3377, respectively. 
The age and metallicity distributions are affected by 



observational incompleteness. To assess the correction that 
likely needs to be applied, we used the functional relations 
between K' , mass, Z and Age using the Bruzual fc Chariot 
( 2003 ) SSP model. Since the HST/ACS and WFPC2 Ri- 
band photometry is 100% complete for all K' detections, 
no additional V, I completeness corrections are necessary. 
We fit power-law and logarithmic functions to the K' vs. 
[Fe/H] and age relations for mass at the 50% completeness 
level at 2 and 14 Gyr and for all metallicities, respectively, 
obtained from the IBruzual & Chariot ( 2003 1 SSP models 



(cf. lines in Fig. [9|. The age and [Fe/H] distributions were 
then convolved with these functional relations to illustrate 
the completeness correction for 50% incompleteness. These 
estimates are shown with light dashed lines in Figure [7] Due 
to the relatively shallow slope of the M — Z relation (cf. 
Fig.[9|, the incompleteness affects equally metal-poor and 
metal-rich GCs, unless if metal-poor GCs tend to be larger 
and less dense which could bias their detection. 

Keeping in mind photometric uncertainties, the age and 
metallicity distributions in Figure[7| seem to indicate the 
presence of multiple GC populations in the sample galax- 



6 COLOUR-MAGNITUDE AND 

MASS-METALLICITY DISTRIBUTIONS 

The colour-magnitude diagrams of GCs in the three merger- 
remnant galaxies is presented in Figure [8] With different 
lines in the left-hand panels of Figure |8] we show [Bruzual &l\ 
Chariot] |2003) SSP model tracks for all metallicities scaled 
to the luminosity of a GC with M = 10^ M© at 14 Gyr, using 
model M / Lv values. The iso-mass lines in the right-hand 
panels in Figure [s] are calculated using the model M/Lv at 
14 Gyr. This comparison shows that the masses of the GCs 
with K' data range from a few IQPMq up to a few 10® Mq. 
The upper and lower solid lines in the left panels labeled 
with 50% completeness indicate the 50% completeness at 
14 and 2 Gyr, respectively, estimated using the SSP model 
prediction for the I — K colour. 

To convert the colour-magnitude diagrams to mass- 
metallicity distributions, we use M/Ly values provided by 
the SSP model to convert from GC luminosities to masses 
using the distance moduli to the galaxies. For comparison 
purposes, we also include info on GCs in the Milky Way 
(hereafter MW GCs). MW GC masses were calculated using 
data from McLaug hlin fc van der Marel| ( |2005[ ) and |Harris| 
(1996). In Figure [9] we show the mass-metallicity distribu- 
tions of all GCs in the sample galaxies. To highlight the 
GCs that are old enough to have developed HBs, we divide 
the GCs in two subsamples: older or younger than 8 Gyr 
(solid blue circles and open red squares, respectively). The 
MW GCs with extendec^HBs, i.e. EHBs, blue HBs, or red 
HBs in the bottom panel are shown with solid (blue) circles, 
(blue) crosses and open (red) squares, respectively. As Fig- 
urel9|shows, all MW GCs more massive than 6 x lO^M© (in- 
dicated with a horizontal dashed line) have hot HBs (EHBs 
or BHBs), except 47Tuc. 

As expected, the GCs in NGC 3610 with spectroscopy 



12 AVhb > 3.5mag;|Lee et al.N2007l 
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Figure 8. Colour-magnitude distribution of GCs (solid circles) in 
NGC 3610, 584 and 3377, from top to bottom, respectively. Large 
circles are GCs with colour error < 0.1 mag. With circles in the 
NGC 3610 panels are indicated GCs with spectroscopic observa- 
tions in the literature jStrader et al.|2003[|2004| . Lines in the left 
panels show 1,2,5 , 10 an d 14Gyr isochrones and isometallicity 
[Bruzual fc Chariot] ( |2003 | tracks for lO^M© and Z indicated in 
the legend. Lines in the right panels mark a constant mass for 
all metallicities at GCs' peak age obtained from the probability 
density distributions in Fig.[7] using the M/Lv-ratios provided 
in the models. The upper and lower solid lines in the left panels 
labeled with 50% completeness indicate the 50% completeness at 
14 and 2Gyr. 



are among the most massive ones [M > 6 x 10^ M©). 
However, among those, only W28 is formally older than 
8Gyr with t = 11.047±i:gti Gyr, and the age of W3 is 
t — 7.9^2 2 Gyr, equal to 8 Gyr within the measurement 
errors. This renders them good candidates for GCs that 
host BHB stars. This further supports the interpretation of 
the difference between photometric and spectroscopic ages 
due to the effect of the BHB in these metal-poor GCs as 
discussed in SectionlOl On the other hand, the younger 
W31 and W6 (t = 63^} g Gyr, [Z/H]= -0.498±g:;i and 
t = 5.9tJjGyr, [Z/H]= -0.610±°oi7t) likely have not yet 
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Figure 9. Mass-metallicity distributions for GCs in NGC 3610, 
NGC 584, NGC 3377 and in the MW (from top to bottom). Solid 
circles and open squares in the sample galaxies panels show GCs 
older and younger than 8 Gyr, respectively. Solid symbols, crosses 
and open squares in the Galactic GCs panel separates GCs with 
very extended, blue and red horizontal branches. Dashed hori- 
zontal line marks the 6 X 10^ Mq, above which all Galactic GCs 
except 47Tuc have BHBs. Labeled solid lines indicate the 50% 
completeness level at 2 and 14 Gyr. 
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developed BHBs, which would explain the consistent photo- 
metric and spectroscopic ages seen in Figure|4] 



7 COMPARING PROPERTIES OF GLOBULAR 
CLUSTERS WITH THOSE OF THEIR HOST 
GALAXIES 

Our aim is to assess information on the star formation histo- 
ries (SFHs) of the galaxies in our sample as determined from 
the properties of their GCs, and compare it to known pop- 
ulation properties of the host galaxies themselves. This will 
not only test for consistency between the two approaches, 
it will also highlight the particular insights gained from the 
combination of deep optical and near-IR photometry of GCs 
in nearby galaxies. 



7.1 Ages and metallicities 

Under the assumption that massive GCs were formed dur- 
ing major star formation epochs in the galaxies' assem- 
bly history, one would expect a relation between the peak 
ages (and metallicities) for the GCs and the (luminosity- 
weighted) ages and metallicities of their host galaxies. To 
test for such a relation, we collected recent age and metal- 
licity measurements for the integrated light of the sample 
galaxies from the literature. These are listed in Table[6l and 
were derived from B-V a.ndB-K colours ( |Li et al.|2007l) 



as well as Lick indices from long-slit spectroscopy by Howell 



1 2005 1 and|Sanchez-Blazquez et al. ( 2006J using the Thomas 
et al. (2003 ) and Vazdekis/MILES SSP models, respectively 



Sanchez-Blazquez et al. ( 2006 1 provide ages and metallici- 
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ties from four different index-index diagrams ([MgbFe]-H;3, 
Fe4383-H^, Mgb-H^ and CN2-H^) and the average of their 
"best nine spectral synthesis model fits" (giving the oldest 
age) . We selected the average value of all these five measure- 
ments, adding errors in quadrature. 

For NGC 3377 there are also age and metallicity esti- 
mates from the SAURON project ( [Kuntschner et al. 2010[) 
at Re, which are 8. it" * Gyr and [Z/H]= -0.23 ± 0.04dex. 
These were derived from comparing Mgb, Fe50, Fe52, and 
Ufj Lick indices to the Vazdekis/MILES SSP models. These 
values are similar to the results of [Sanchez-Blazquez et aL] 
1 2006] ) for NGC 3377. The SAURON values are not included 
in Table[6] because the other two galaxies were not in the 
SAURON survey. An interesting result on NGC 3377 from 
the latter survey is that its inner 20arcsec (~ 0.7 kpc) re- 
gion features a relatively young stellar population with high 
metallicity located in a rotating structure with disc-like kine- 
matics ( Kuntschner et ar]|2010 l. This further supports the 
accretion/interaction-driven formation history of NGC 3377. 

Probing for a correlation between the luminosity- 
weighted stellar age and metallicity of the galaxy and its 
GCs, we plot the peak age and metallicity of the GCs versus 
host galaxy age and metallicity in Figure fTOl It is clear that 
the peak age of the red (i.e. metal-rich) GCs in all galaxies 
is indeed correlated with the luminosity-weighted mean ages 
and metallicity of the host galaxy. 

The upper panel of Figure [To] suggests that the 
(luminosity-weighted) ages of the sample galaxies are typi- 
cally younger than the mean ages of their bright metal-rich 
GCs. While this difference is likely due in part to systematic 



Figure 10. Correlation between ages (top panel) and metallici- 
ties (bottom panel) of red globular clusters and integrated light 
of their host galaxies. Different symbol types represent different 
literature sources for the (luminosity-weighted) mean age and 
metallicity of the galaxy, as indicated in the legend. Dashed line 
is the one-to-one relation. 



effects related to the use of different SSP models, this sug- 
gests, at face value, that the galaxies had a prolonged star 
formation history since the major epoch(s) of massive GC 
formation. One scenario that could cause this effect is that 
the star formation occurring more recently than the era in- 
dicated by the peak age of the metal-rich GCs featured lower 
star formation rates. This would form fewer high-mass GCs 
and thus cause fewer GC detections at a given detection 
threshold, while still causing higher surface brightness for 
the integrated galaxy light, thus yielding younger ages for 
the latter. This scenario can be tested directly by probing 
less massive red GCs with deeper near-IR imaging data in 
the future, e.g. using the NIRCAM instrument on the James 
Webb Space Telescope. 



7.2 Reconstructing galaxy SFH through its most 
massive GC 

Following recent development on the relation between 
galaxy star formation rate (SFR) and the luminosity of its 
most massive cluster at a given epoch, we attempt to derive 
the galaxy star formation history (SFH) in this section. 

The formation efficiency of young massive clusters that 
can form in a galaxy has been shown to scale with the galaxy 
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Table 6. Literature values of the age and metallicity of the galaxies in our sample derived from B — V,B — K colours ||Li et al.|2007| 
(columns 2 and 3), and from integrated light spectroscopy based on Lick indices by |HowellH 2005 > based on [Thomas et al.| | |2003| l SSPs 
(columns 4 and 5) and by [Sanchez-Blazquez et al.| f2006) using Vazdekis/MILES SSPs (columns 6 and 7). In columns 8 and 9 we show 
the peak age and metallicity of the metal-rich GC subpopulation. 



Age [Gyr] 



[Z/H] 



Age 



[Z/H] 



Age 



from (Li et al.l2007l 



(2) 



(37 



I Howelll2005 I 



(41 



15) 



I Sanchez-Blazquez et al. 



[Z/Ifl I GCs Age [Z/H] 

from this study 
IT) (8) (9) 



1.60 ±0.22 0.389 ±0.034 1.7 ±0.1 0.76 ±0.16 

3.13 ±0.72 0.279 ±0.060 2.4 ± 0.3 0.61 ± 0.04 6.2 ± 0.69 0.125 ±0.031 
4.08 ± 1.00 0.061 ±0.050 3.5 ±0.8 0.30 ± 0.06 8.9 ±1.58 -0.010 ± 0.068 



4± 1 
6± 1 
9± 1 



0.02 ±0.15 
0.19 ±0.13 
-0.28 ± 0.16 



SFR ( [Larsen fc Richtler|[2000|. Since the SFR is higher for 
higher gas densities ( Kennicuttj 1998 1 , the formation of mas- 
sive clusters becomes more likely as well (Larsen & Richtler 
|2000[ [Ehnegreen fc Efremov||1997[ ). However, the formation 
of a massive cluster can also be subject to statistical fluctu- 
ation resulting from a size-of-sample effect in which higher 
SFR galaxies will form more GCs and thus sample better 



the high-mass end of the cluster mass function ( Billett et al 



2002 Larsen 2002 1. Under the assumption that the observed 



relation between the brightest cluster and the galaxy SFR 



( Larsen 2002 1 is explained by a physical limit for the pos- 



sible range of cluster masses regulated by the galaxy SFR, 
[Weidner et al.| ( |2004[ ) derived relations between the galaxy 
SFR and the initial (embedded) mass of its most luminous 
cluster, which allows one to recover the galaxy SFR. To ob- 
tain correctly the galaxy SFH, the present-day cluster mass 
has to be corrected for stellar evolution and dynamical mass- 
loss to obtain the mass of the cluster at the time of its for- 
mation ( jMaschberger fc Kroupa 2007). However, the most 
luminous cluster in a galaxy is not always the most mas- 
sive one. For starburst galaxies, the most massive cluster 
is typically young as shown in Monte Carlo simulations by 
Bastian (20081. Conversely, if the duration of the Monte 



Carlo experiment is extended to few Gyr, the oldest clusters 
are usually the brightest ( |Gieles|["2009| ). Nevertheless, |Bas- 



tian 



( 2008 1 concludes from his Monte Carlo simulations that 



the observed relation between galaxy SFR and its brightest 
cluster accurately reflects the recent SFR. Thus, Bastian 



employs the M, 



brightest 
V 



VS. SFR relation from 



Weidner et al. 



( 2004 1 and corrects for stellar evolution, but not for dynam- 
ical mass loss (expected to be small for massive clusters, 
since the tidal/dynamical mass-loss is inversely proportional 
to cluster mass; Baumgardt fc Makino|2003 1, the masses of 
the brightest GCs associated with a given epoch for a sample 
of major merger galaxies to derive their SFHs. Among the 
sample of their galaxies is NGC 3610 for which they used 
the Mv magnitudes from Whitmore et al. ( 2002 1 and de- 



rived a peak SFR of 218 Mq yr"\ Here, we follow the same 
approach to assess the SFH of the galaxies in our sample as 
follows. 

Using the photometric masses of our clusters, we select 
the five most massive GC in red three age bins, as shown 
by vertical dashed lines and indicated with arrows in Fig- 
ure |ll[ As several studies have shown that metal-poor GCs 
and metal-rich GCs trace different components of early-type 
galaxies, i.e., the halo and bulge/spheroidal components, re- 
spectively (e.g., Kissler-Patig et al. 1997 Kundu fc Whit- 
more |p98l [Goudfrooij et al.||2007| 7 we distinguish between 
the metal-poor and metal-rich GCs using I ~ K' = 1.4 as 



a dividing colour as before (cf. Section |5.3[ ). The mass of 
the most massive GCs in those age bins is then used to 
calculate the galaxy SFR using Eq. 7 of [Weidner et al.| 
(2004). To correct the mass of the GCs for mass- loss due 
to stellar evolution, we used predictions from the [Bruzual] 
fc Chariot] ( [2003 j ) SSP models. We used the SSP M/Lv at 
lOMyr to estimate the clusters' initial masses, as in [Bas-[ 
(|2008|). In Figure[l2| we show the derived galaxy SFRs 



tian 



from the five most massive GC in three age bins (t < 4 Gyr 
<t<8Gyr, and t>8Gyr), where their age positions in the 
figure are centered around the midpoint of the respective age 
bin. The peak SFR we derive for NGC 3610 (234 Mq/ yr) i s 
consistent with the 218M0/yr derived by [Bastian] ( |2008| ), 
taking the difference in assumed distance moduli into ac- 



count they adopt to NGC 3610 from |Whitmore et al.| ( |2002 l 
(A(m — M) = 0.32 mag), as well as the adopted by them 
solar metallicity for the most massive GC. For the three 
sample galaxies, Figure[T2] shows that the peak star forma- 
tion was reached at the epoch of formation of the oldest, 
metal-poor GCs. Since then, dissipative, gas-rich episodic 
(merger/accretion) events triggered the formation of the 
younger and more metal-rich stellar populations. 

Similarly, we calculate the SFRs of the Milky Way and 
M31 from their five most massive old (t > 8 Gyr) GCs (cf 
Figs. [5] and 111, which are representative for the SFRs at 
the time of assembly of their spheroid components (bulge 
and halo). We believe this yields a useful depiction of the 
uncertainties involved with the use of the most massive GC 
as a probe of the SFRs within these age bins. We do not 
expect observational bias effects in the SFR-values due to 
incompleteness (all measured "brightest" GCs are above the 
90% completeness level, see Fig. 11 1, however we can not ex- 



clude the possibility of missing a massive GC which might 
be located outside of the NIRI and HST WFPC2 or ACS 
fields (cf. Fig.[TJ . Effects on the SFRs due to statistical fluc- 
tuations in GC mass distribution ( [Ma schberger fc Kroupa 



[2007| [Kruijssen fc Cooper|20lT| ) must be also kept in mind. 

The result of our analysis supports the application of 
the most massive GCs being a good tracer of the galaxy SFH. 



8 CONCLUSIONS 

We use new Gemini/NIRI K'-ha,nd imaging in conjunction 
with existing optical V, I HST/WFPC2 and ACS photom- 
etry of bright GCs in NGC 3610, NGC 584 and NGC 3377 
three early-type galaxies that feature signs of recent galaxy 
interactions. The addition of the I — K colour effectively 
breaks the age-metallicity degeneracy present in the V — I 
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Figure 11. Mass as a function of age for the GCs in NGC 3610, 
584 and 3377 from top to bottom. Solid and open symbols indicate 
clusters bluer or redder than I — K = 1.4mag, i.e. metal-poor or 
metal-rich. Vertical dashed lines show the age bins within which 
the most massive GCs were selected to calculate the galaxy SFRs 
shown in Fig. |12| Dashed curves show the 50% and 90% complete- 
ness levels estimated from the SSP models for all metallicities. 
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colour. We interpolate between Bruzual & Chariot ( 2003 \ 



SSP model tracks to derive photometric ages, metallici- 
ties and masses for all GCs in the sample galaxies. For 
the massive GCs {M ^ 6 x 10^ Mq) in NGC 3610 with 
available spectroscopic Lick indices measurements from the 
literature JStrader et al. 1 120031 [20041, we find that the 



photometric ages are older by ~ 2 Gyr. This age differ- 
ence is larger for the metal-poor GCs, while the photomet- 
ric and spectroscopic metallicities are in excellent agreement 



(Sect. 5.1 1. We argue that this is most likely due to the pres- 



ence of hot blue HB stars: A comparison with HB prop- 
erties of Galactic GCs shows that all Galactic GCs with 
M^6x 1O^M0 have blue HBs, except the metal-rich GC 
47Tuc. This suggests that all NGC 3610 GCs with spec- 
troscopy found to be older than ~ 8 Gyr are most likely to 
possess a blue HB as well as the majority of metal-poor GCs 
with 7Vf>6xl0'^ Mq (in any galaxy) . 

We show that the peak value of the age and metallic- 
ity distributions of the sample galaxies GCs is correlated 
with the luminosity-weighted mean age and metallicity of 
the host galaxy. The GCs' age and Z distributions are broad 



(Sect. 5.31, indicating a prolonged star and cluster forma- 



tion history. We reconstruct the galaxies' SFHs from the 
GC age distributions, using the recently observed correla- 
tion between the luminosity of the most massive GC at a 



given age and the galaxy SFR (Sect. 7.21. 

Our results support a scenario in which the star for- 
mation rate of the three candidate intermediate-age galax- 
ies in our sample peaked at the epoch at which the oldest 
GCs (i.e., the metal-poor ones) formed. Subsequent dissi- 
pative events (i.e., interactions involving gas-rich galaxies) 
led to the formation of more metal-rich GCs up to a few 
Gyr ago. The SFHs of our sample galaxies as determined 
from their GCs suggests that the SFR was highest at the 
earliest times (> 10 Gyr ago in all three cases. The latter 
is in large part due to the presence of massive metal-poor 



Galaxies' SFHs from V, I, K' photometry of Globular Glusters 17 



GCs which are likely associated with the early build-up of 
the galaxy halos. The strengths of subsequent star formation 
events occurring over the last ~ 10 Gyr varied significantly 
from one galaxy to another. Among the three galaxies in our 
sample, NGC 3610 showed the strongest SFR in all age bins, 
especially at the youngest age bin (<4 Gyr ago). This is con- 
sistent with NGC 3610 being the most luminous galaxy in 
our sample as well as having the youngest spectroscopically 
determined age. 

The ability to obtain such direct insights on the SFHs of 
galaxies renders the study of galactic GC systems in the op- 
tical vs. near-IR a promising tool to trace galaxy formation, 
merging and star formation history. 
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